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SUMMARY 

In this paper, maximum likelihood estimation for distributed parameter models of large 
flexible structures has been formulated. Distributed parameter models involve far fewer unknown 
parameters than independent modal characteristics or finite element models. The closed- 
torm solutions for the partial differential equations with corresponding boundary conditions 
have been derived. The closed-form expressions of the sensitivity functions lead to highly 
efficient algorithms for analyzing ground or on-orbit test results. For illustration of this 
approach, experimental data of the NASA Mini-MAST trust have been used. The estimations 
of modal properties involve its lateral bending modes and torsional modes. The results show 
that distributed parameter models are promising in the parameter estimation of large flexible 
structures. 
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INTRODUCTION 

Large spacecraft structures, such as, Solar Array Flight Experiment, Mim-MAST (CSI) 
employ large complex trusses in their constructions. Because of their large size and because of 
gravitational loads, it is not possible to determine with suitable accuracy their structural dynamic 
characteristics from ground-based testing of a full-scale prototype. Analysis of on-orbit response 
data will be necessary. In recent years, numerous variations of system identification methods 
have been developed. Unfortunately, current approaches to parameter estimation cannot handle 
the complex models foreseen to be necessary. As the number of modes increases, the accuracy 
will be decreased and the complexity significantly increases. 

Two distinct approaches to the solution of large space structural parameter identification 
problems have emerged. One is the lumped parameter approach, the other is the distributed 
parameter approach [1-5]. The obvious fact is that by far the most effort put into the study of 
the identification problem has been based on the lumped parameter model. With the increase 
of the number of modes, the number of parameters increases rapidly for the lumped parameter 
approach. However, for the distributed parameter approach, the total number of parameters 
needed to be identified almost keeps the same. 

Distributed parameter model is based on the partial differential equation (PDE). In this 
approach, instead of identifying the modal frequency, damping and mode shape deflection of 
each mode, only the coefficients of the partial differential equation and initial conditions need 

to be estimated. 

This paper intends to create a very simple distributed parameter model, combining the use of 
maximum likelihood estimation technique (MLE), to identify the modal characteristics of NASA 
Mini-MAST truss which is treated as a cantilevered beam with two concentrated masses. Wave 
equation and Bernoulli-Euler beam equation are introduced to describe the torsional and bending 
behavior of the Mini-MAST truss. Proportional damping will be taken into account simply by 
adding a damping term in the PDE’s, which is proportional to velocity (6]. 
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A series of closed-form solutions of the PDE’s will be used to match the measurements of the 
vibrations. Based on the optimal fitness between the measurement response and the theoretical 
response in the sense of maximization of likelihood function, the coefficients in the PDE’s and the 
parameters dependent on the initial conditions are estimated. The closed-form expressions of the 
sensitivity functions are derived, which carry out the innovations of the unknown parameters in 
the iteration process. The comparison of the results with other methods shows that the proposed 
approach is promising in the parameter estimation of large flexible structures. 

NASA MINI-MAST (CSI) TESTBED 

NASA Mini-MAST (see Fig. 1) is a generic space truss built primarily for research in the 
areas of structural analysis and vibrational control [7]. Mini-MAST is deployed vertically inside 
a high-bay tower, cantilevered from its base on a rigid foundation. The total height of the truss 
is 20.16 meters, containing 18 bays (1.12 meters each) in a single-laced configuration with every 
other bay repeating. The deployable/retractable truss beam has a three- longeron construction 
forming a triangular cross-section with points inscribed by a circle of diameter 1.4 meters. The 
beam has three member types: longerons, battens, and diagonals. Longerons are parallel to the 
beam axis and provide beam stiffness and strength in bending. Battens are in the beam face 
planes and provide beam stability. Finally, diagonals, also in the beam face planes, provide beam 
stiffness and strength in torsion and shear. Two instrumentation platforms, holding actuators and 
sensors for CSI control experiment, have been installed at Bay 10 and Bay 18 (beam tip). These 
additional components have a significant effect on the structural dynamic characteristics. 

In this paper, the real Mini-MAST truss is treated as a cantilevered beam with two 
concentrated masses at Bay 10 and Bay 18 respectively (Fig. 2). The continuum model is 
equivalent to the real truss in the sense that both have the same dynamic properties, say, natural 
frequency, damping ratio and mode shape. To keep the equivalency, the structural parameters 
of the equivalent must be set up properly. AH these parameters will be divided into two types. 
Some physical quantities of the real structure, such as length, weight, etc., are called unadjustable 
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parameters which are assumed to be known a priori. Another type of parameters, such as the 
composed parameters appearing in the PDE’s, will be adjusted in the estimation process. This 
type of parameters is called adjustable parameters, which will be the elements of the unknown 
parameter vector in the maximum likelihood estimator. The parameters of the equivalent of the 
Mini-MAST truss is listed in Table 1. 

The control inputs of the system are three orthogonal torque wheel actuators located at the 
top platform (Bay 18). Output pulse responses were obtained by applying single pulse signals 
at each input channel. Twenty seconds of output pulse response were collected for each input 
channel. A sampling frequency of 50 Hz is used. Two sets of data [8] are selected for our 
analysis here. The acquisitions of the selected data are recounted as follows. The first set 
of data was obtained from the measurement of the rotation rate about the x-axis, which was 
measured by one rate gyro mounted at the tip platform. The second one was obtained from 
the measurement of y-direction displacement, which was measured by one displacement sensor 
installed at Bay 18, mounted parallel to the fiat face on the comer joints of the structure and 
positioned to measure deflections normal to the face. The locations of the actuators and sensors 
concerned are shown in Fig. 3. 


MODELLING 


In this paper, the analysis of modal characteristics of NASA Mini-MAST truss involves 
torsional and bending motions. The damped torsional vibration is described by the wave 


equation 


o-o oe , o'-o 

+ c ~5i ~ k d? 


(i) 


where 9{x,t) angle of twist, Jt, = plb moment of inertia of the beam, k = Gib torsional rigidity, and 
C damping constant of proportionality. Two parameters a and b, which relate the coefficients 
of the PDE, are defined by 


a 


2 



( 2 ) 
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The lateral bending vibration is described by the Bernoulli-Euler beam equation with proportional 
damping term 


d 2 y 9>J , d*y 

m dt* +c dt + k dx i 


- 0 


(3) 


where y(x,t ) lateral displacement, m = pA mass per length of the beam, and k = El lateral 
bending rigidity. Two parameters a and b are also defined by 


Hie solutions of Eqs. (1) and (3) can be written in a generic form (see Appendix) 


V (M) = r ' ( x ) e ( A * cos “dJ + D * sin u d ,t) (5) 

I 

where, y(x,t) represents either angle of twist (0) or lateral displacement y, A x and B t the 
coefficients dependent on the initial conditions, u) Ht and yj 1 — the natural frequency 

and damped natural frecj uency , respectively, and £/ is the damping ratio. All these modal 
properties are related to the parameters a and b by 


w «. = “ k l u d, = \J {nklf - b 2 , e, = (6) 

where, q represents the order of power, with q — 1 for torsion and q — 2 for bending. From Eq. 
(6), the solution Eq. (5) can be expressed in terms of the parameters a and b, 


V ( x , t) — y y Vi (x)e hl .4, cos t \J (nA.?) 2 -b 2 + B, sin t yj (a A.?) 2 - b 2 


(7) 


l,(x) in the solutions (5) or (7) are the eigenfunctions given as follows. For torsional 
equation: 

li(*)= l 


b( W + 1 


sin k{X 


{ + sin k ' x 


0 < x < X|q 
X[(j < x < l 


( 8 ) 
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where. 


»(M) = 


sin W + +M COS A;,/ 


cos A-,/ - -jjT',/ siu k x l 

For Bernoulli-Euler equation 

n (sh A,x - sin ki x) + (cli k,x - cos k, x) 


Yi (*) 


0 < x < xio 


a, [/i , (k,l) sin k,x \ h.\(k,l ) cos k t x + shk t x] p 
fj , [/13 ( k,l) sin k,x + h‘i (A,/) cos k t x + chk t x] 


arm < X < l 


where. 


and 


( 9 ) 


(1,0 = „ (fcMfcw-w ‘W mv> = «<w) - 

l„(k,l) = 93(M)+2jftWc«,kil cftl-,1; ft, (1.0 =9t(l.O - Mil 

m (til) = os ft, I <./<ft,l+ sin 1,1 si, ft, I, n (ft.o = cos ft, I Cftft.l - sin 1,1 sftft.l 

, (; ,(|,1) ... ,-„»ft,l si, ft, I I sin ft, I , 1,1,1. 9, (ft.O - cos ft, I sftft.l • siiift.l cftft.l 
and a;, /7, and r, are the modal participant coefficients. 

Finally, the most important quantities are the eigenvalues Ki s of the system. They are 

roots of the corresponding characteristic equations given as follows, 

For torsional equation: 

kl sin 2 Ax io ■+ — j^-kl g(kl) sin kxm — 1 ^ 

5 . 7 Jb 


For Bernoulli-Euler equation: 



Cli X 

-SHX -CHX 

sax 

n ] SNx^h i csx 

h$S N X ■+ h 2 C ^ X 

SNX CSX 

hi CSX -In SNX+CH X 

h j C.’ 5 X' - hj S N X + S11X 

CSX -CHX -SNX-SHX 

iCH X 

l SHX 

-it •// x + ^ SN X - SH X ) -a sm x + ^ *•'( 


(ID 


where, 


SNX = sin kx\u,CSX = cos Ax 10) SHX = sh Ax 10) and CHX - ch Ax 10 
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Maximum Likelihood Estimator 

Assume that the outcome Y of an experiment depends on an unknown parameter vector 9. 
We want to estimate the best value of 9 according to the observation Y. One of the advanced 
technique is based on the Maximum Likelihood Estimate (MLE) principle. The essence of 
the MLE method is to maximize the conditional probability density function, i.e., so-called 
likelihood function, P(Y\9). 

Suppose we have the measurement response sequences t/(l), y(2 ) , y (in). The matrix 
Y m consisting of all measured outputs is introduced, Y,„ — [y(l) , y (2) , y (m)]. If the 
probability distribution of Y, n has a density P[Y m \9], it then follows from the definition of 
conditional probability that 

771 

P [Y m \6] = P [y (1) , y (2) , .... y ( m ) \9\ = JJ P [y (i) \ Y,„ u 9} (12) 

1 = 1 

If the assumption of Gaussian distribution is taken, the likelihood function has the fonn of 

m i f i 

m = I’tn.M - n - >.*> -sW - M» 

■ R~ l Mil* - M) - !/(■!' - M)] j (13) 

where y u (i|i — 1, 9) is the norminal response calculated by using 9 0 . If the constants are ignored, 
we have the log-likelihood function 

J(9 ) = - hi L(9) 

m 

= 2 _ M) - y(*l* - M)] ; ^ _1 (y (*1* - M) - y(*l* - M)] + 

Z i=l 

(14) 

Linearizing y(i\i —1,9) with respect to the unknown parameter vector 9, we have 

y(*|*- M) = j/«(*|t- M) + (V0?/,)(0 - 9 0 ) (15) 
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where V*y, gradient of y with respect to 0, and 0 O is nominal 0 vector. Substituting Eq. (15) 
into Eq. (14) for J(0) and setting - 0, we may obtain. 

“I 

\T _ 'Ti \\ (16) 


8 = 8 a + 


YAVm) R l (^eVi) 


i=i 


5Z(V^yi) ii -1 (Vi — i/i«) 

L.=i 


In this paper, the unknown parameter vector will be 

0 = [a, b , A U A 2) ...,A n , B lt B 2 , 5./ 


The relationship between the modal properties and the unknowns is given in Eq. (6). Thus 
we can obtain the modal properties through the solution of the PDE as long as these unknown 
parameters are determined. The gradients of y(*,l) with respect to the unknowns in Eq. (16) 
can be derived from Eq. (7) simply by taking derivatives. 


2a. ( X t) = 'VYj(x)e- it - piL = - [-^sin</(qfcf) 2 -6 2 + B t cos b 1 

Oa ’ V JtnLrf-V L 




2a t)^yv, (0 f<T 6< ( -Ai cost^afc?) 2 -"^ - D, sin t/Kf " 

- B, costs/ (aktf - V | 08) 


Vw - ^ 


_5y_ 

5.4, 


9B, 


(x-,0 = YY*)e' bt cos t\J{aktf ~b 2 


(x,t) = Y,{x)e- bt xmtyf (ak v t ) 2 - b 1 


(19) 


(20) 
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Analysis of Modal Properties 

As mentioned before, the procedure by using distributed parameter approach to analyze the 
modal properties is quite different from that of the lumped parameter approach. First, based on 
the unadjustable structural parameters, the eigenvalues k, and the eigenfunctions K, (x) (mode 
shapes) of the system can be determined through the solutions of the corresponding characteristic 
equations before the estimation iteration starts. Second, only the coefficients of the PDE’s and 
the parameters relevant to the initial conditions need to be estimated rather than the modal 
frequency, damping and mode shape deflection of each mode. These two characteristics greatly 
decrease the number of unknown parameters and speed up the iterative process. 

Solution of the PDE is in the form of infinite series mathematically, so no modal truncation 
problem is involved theoretically. However the contributions to the response are always so small 
for the higher frequency modes that only the first several modes (five in this paper) are used 
in the analysis. 

It is noted that the modal coupling must be considered. Because of the eccentric properties of 
both the tip-mass and the tip-actuating pulses, the lateral bending vibration will be excited while 
torsional vibration exists and vice versa. However, the experimental data show that bending 
modes are hardly recognized in the torsional measurement, and the first torsional mode appears 
clearly in the bending measurement. So the first torsional mode is included in the analysis of 
the bending vibration. 

From the experiment data it is hard to get any a priori information about the initial conditions 
of the response, that is, the initial values of .4, and D x for iteration. Fortunately the initial values 
of .4, and B % do not affect the convergency significantly, so they are chosen arbitrarily. However, 
the initial values of parameters a and b are very important to the convergency. From the results 
of finite element analysis, we can determine the initial values of equivalent stiffness first, then 
proceed to reckon the initial values of the parameters a and b. 
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Iterative accuracy is controlled by the innovation residual of the unknown parameter vector 
0, which is defined by 

x 1/2 

e# = 

where, p the number of the unknowns, 0(i) the ilh element of the unknown parameter vector 
0, and n is the number of successive iteration. In the algorithm, eg < 1(T 7 is referred to as 
the criterion which controls the iteration. 

Table 2 shows the comparison of the estimated frequencies obtained from Finite Element 
Analysis (FEA), Eigensystem Realization Algorithm (ERA) [9] and Distributed Parameter Algo- 
rithm (DPA). Most modes are comparable to each other. The fifth bending mode from DPA is 
extremely higher than that of the other approaches. This is due to inadequacy of the Bernoulli- 
Euler beam model used. Because the rotary inertia and shear deformation of the beam are 
neglected, the Bernoulli-Euler beam model produces much higher frequencies in the high fre- 
quency range. In order to improve the accuracy for high frequency, Timoshenko beam model 
is proposed for further investigation. 

Figure 4 shows that the reconstructed responses obtained from the estimated parameters and 
the measured responses have a reasonably good fitness. 

CONCLUDING REMARKS 

This paper proposes a distributed parameter model for the analysis of the modal 
teristics of NASA Mini-MAST truss. Wave equation and Bernoulli-Euler equation have been 
used to describe the torsional and bending vibrations respectively. Closed-form solutions of the 
PDF's are derived. By using the Maximum Likelihood Estimation method to provide the opti- 
mal match between the experimental data and estimated responses, the coefficients in the PDF’s 
and the parameters dependent on the initial conditions are estimated and the modal properties 
can be further determined. The results are comparable to those from other approaches. The 
estimates of bending modes in the higher frequency range is expected to be improved by using 
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Timoshenko beam model. Because the number of unknown parameters is greatly reduced in 
the distributed parameter model and the maximum likelihood estimation is feasible based on the 
derived closed form solutions of the PDE’s, the proposed approach is particularly attractive for 
its less computational burden for the large flexible structures. 
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FIGURE LEGEND 


Fig. 1 Sketch of NASA Mini-MAST Truss 

Fig. 2 Equivalent Cantilevered Beam 

Fig. 3 Mini-MAST Sensor and Actuator Locations 

Fig. 4 Comparison of Reconstructed and Measured Responses 
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Table 1 Structural Parameters of Equivalent Continuum of NASA Mini-MAST Truss 

A. Unadjustable Parameters 

Length of the Beam L = 66.24 ft. 

X-Coordinate of Bay 10 x 10 = 36.80 ft. 

Radius of Gyration of the Section r = 1.6237 ft. 

Mass per Length pA = 0.1076 slug/ft. 

Ratios of Weights: 

Bay 10-body/Beam: = 0.4760 

Tip-body/Beam: = 1.4547 

Ratios of Moments of Inertia: 

Bay 10-body/Beam: = 0.6206 

Tip-body/Beam: ^ = 0.6206 

B. Adjustable Parameters: Initial Values for Iteration 

Longitudinal Stiffness EA = 10,530,000 Lb 

Bending Rigidity El = 27,760,000 Lb.ft. 2 

Torsional Rigidity GI p = 1,970,000 Lb.ft. 2 


<994 
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Sketch of NASA Mini-MAST Truss Fig. 2 Equivalent Cantilevered Beam 













Response (ft.) Response (ft.) 



Fig. 4 Comparison of Reconstructed and Measured Responses 
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appendix a 

Solution to the Bernoulli-Euler Beam Equation 


In using distributed parameter approach to identify the structural systems, one of the most 
important procedures is to solve the partial differential equation. 

Here, only solution to the Bernoulli-Euler beam equation will be provided. Solution to the 
torsional vibration equation can be derived by following similar procedures. 

The Bernoulli-Euler beam equation describing lateral bending is 

d*y 1 d 2 y 


= 0 


dx* a 2 dt“ 

where, a 2 = fj. The general solution to Eq. (A.l) may be expressed as 

y (x, *) = Y i ( x ) [* 4 > cos + B i s * llu ’>^] 


(A.l) 


(A. 2) 


(A. 3) 


where, Y t (x) are the eigenfunctions which are of the form 

Y, (x) = Ci sin fcjz + C> cos k t x + ( $ s h k,x t C\ch k,x 

Now we derive the specific form of the solution for the equivalent Mini-MAST truss — 
cantilevered beam with two lumped masses. The procedure consists of three steps as follows. 

Step 1 For the right segment, i.e. x = xio ~ l- 

The PDE for the lateral vibration of the right segment is 

Art , J_ Art = n (A. 4) 

dx 4 ' a 2 dt 2 

The boundary conditions for the free end (z = l) are: 

'^(M) = o (A5) 


After the separation of variables 


yn(x,t) = Y R (x)T(t) 


A.l 
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we have a \' H - ODE, 


YH"(x)-k*Y R (z) = 0 

with B.C.’s 

I VI (0 = 0 

The general solution to Eq. (A. 6) and the corresponding derivatives are as follows, 

} (a*) = C\ sin kx + C 2 cos kx + C$shkx C^chkx 

1 ^ (a) — A ( C j cos kx — C> sin kx -f C^chkx f C^shkx) 

^ it ( ,f ) ^ (~^'l ~ ^2 cos Aia -f C^shkx C\chkx ) 

fc ^7” ( u ‘) = ^ 3 (“C'l cos At.t + C 2 sin kx + C^chkx + C^shkx) 

By using the B.C., we have 

f C i = Cj ^f/i (kl) + -W^kl cos klshkl^ f C4 ^3 (A:/) + 2-j^A:/ cos klchkl^j 

\ C, = C-J (</ 4 (A,/) - 2{^A/sinA/sAA/j + C* 4 (</ 2 (A/) - 2]fyklsin klchkl) 

where, 

' </i (A7) = cos klchkl + sin klshkl 

02 (H) — cos klchkl — sinA:/s/iA,7 
i 

( J:\ (kl) = cos klshkl -f sin klchkl 
, g A (kl) = cos klshkl — sin klchkl 

Thus, F /f (x) and its derivatives can be expressed in terms of the coefficients C 3 and C 
^ /< ( ;i ‘) = C3 [/r 1 (A7) sin Ax + A 4 (A7) cos kx + sAAx] 

-(- C 4 [A3 (A7) sin Ax -| A 2 (A7) cos Ax f cAAx] 

^ It (-r) — C3 A [A] (A7) cos Ax — A 4 (A7) sin Ax + cAAx) 

I C, K (Aj ( kl) cos kx - h-> ( kl) sin Ax + sA Ax] 

]}i ( — C'sK 1 [ — A j (A/Jsin Ax - A 4 (kl) cos kx + sAAx] 

+Cih 2 [—A3 (A/) sin Ax — A 2 (A7) cosAx -f cAAx] 

17 ( x ) ~ C3/1 3 [— A j (A7) cosAx + A 4 (A/) sin Ax + cAAx] 

1 +C 4 /v 3 [— Aj (kl) cos kx -f A 2 (A7)sin Ax -f sAAx] 

900 a.2 


(A. 6) 


(A. 7) 


(A. 8) 


(A. 9) 


(A. 10) 


(A.ll) 



where, 


/( l ( kl ) = (j\ (kl) + 2|^ kl cos klshkl 
/*2 (kl) — (Ji (kl) — Zipfckl sin klchkl 
hs (kl) = {/i (kl) + 2 yfckl cos klchkl 
k /i 4 (fcl) = g 4 (kl) - 2%fckls\n klshkl 


Step 2 For the left segment, i.e. x = 0 ~ xio- 

The PDE for the lateral vibration of the left segment is 

O'VL , 1 &VL 


Ox 4 + a 2 


= 0 


The boundary conditions for the fixed end (x = 0) are: 

Jvl(0,<) = 0 

l% t (0,i) = 0 


After the separation of variables 


VL(x,t) = Y L (x)T(t) 


we have a Y), — ODE 

YZ"(x)-K*Y L (x) = 0 


(A. 13) 


(A. 14) 


(A. 15) 


with B.C.’s 

r n (o) = o 
i n (o) = o 

The general solution to Eq. (A. 15) has the form of 

Y'l (x) = di SU1 kx + d> cos kx + d^shkx \ d 4 chkx 


(A. 16) 


(A. 17) 


By using the B.C., we can express Y L (x) and its derivatives in terms of the coefficients <i 3 


and d 4 

Y ij (x) — d,i (shkx — sin kx) T d 4 (chkx cos kx) 

Y[ (x) — dj K (chkx — cos kx) T d\ A (shkx T sin kx) 

V r " (x) — </( A “ (shkx -f- sin kx) + d\ A ” (chkx I cos kx) 
Y\" (x) = dj K* (chkx + cos kx) + d 4 A' 3 (s/iA:x - sin kx) 


A. 3 
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3 I he compatible conditions at the common point, x — xio, between the two segments. 

I. The lateral displacement must be the same, i.e., V/,(x| 0 ) = Yn(x m). Thus, 

di(shkx\o — sin t rio) T <A» (c/iA.xio — cosAxio) 

= Ci [A»t (A:/) sin kx\ a + h 4 (kl) cos A-xjo + s/iA-xio] 

YC.\ [hi (kl) sin kx 10 + hi (kl) cos ixjo + chkxm] (A. 19) 

II. The slope of the center line of the beam must be the same, i.e, Y[ (x 10 ) = Y' n (x 10 ). Thus 

(chkx io - cos A:x J0 ) + d 4 (shkx iti + sin Ax 10 ) 

= C'i [/tj (kl) cos fcxjo — h 4 (kl) sin kx\o + chkx\o] 

+C4 [hi (kl) cos kx\o — fi 2 (kl) sin A.*X]o -f shkxio] (A. 20) 

III. The bending moment must be the same, i.e, r"(x 10 ) = r"(x 10 ). Thus, 

<ij (shkx\o + sin kxm) T d 4 (chkxio + cos A-xjo) 

— [ h 1 ( kl) sin A*x j q h 4 (kl) cos kx\o T s/iA’xnj] 

+C'.j [ — /*3 (A/) sin A.X10 — I 12 (kl) cos fcxio -f c/tA-X)o] (A. 21) 

IV. Because of the inertia of the lumped weight VV10, the shear has a jump at x = x lu , i.e., 

EIY’i (x I0 )T - EIYft (x 10 ) T + ^Y L (x 10 ) T 

9 

which yields 

Y'i" (*■<>) = Y'n (*tu) - ~-lk A Y L (x 10 ) 

*'b 
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Thus 


c h (chkx 10 + cos kx ]0 ) + </ 4 (shkx \ 0 ~ sin fcxio) 

= Ci [— /i| (kl) cos^xiu + hi (kl) sin fcxio + c/ifcxio] 


+Ci [—hs (kl) coskx\o + /i 2 (kl) sin fcxio + s/tfcxjo] 


Ww 

W b 


kl [ds (shkx 10 — sin fcxio) + <h (chkx io — cos fczio)] 


(A. 22) 


Then a set of equations forms from Eq. (A. 19) to Eq. (A.22) 


r C’ishkx io I Cichkx [o — dishkx\ii d\chkx io 0 

C';i (h i sin kx io -(- /i 4 cos Aixjo) + C\ (hi sin kx 10 + h > cos kx 10) 


< 


T</j sin kx m d\ cos kx jo = 0 
Cs (h\ cos kx id — hi sin kx m I chkx\o) + Ci(hi cos kx 10 

— h> sin kxm + shkxm) + 1/3 (cos A:xio — chkx iu) 

- d\ (sin kx m T shkx 10) — 0 


Ci( 2 chkxio) + Ci ( 2 shkx\u) + d s (- 2 chkxm) + Tfiffkl (sin fcxiu shkx io)j 


+d 4 (— 2 shkxio) + j^)kl(< 


= 0 




(A. 23) 


The condition for set (A.23) having non-trivial solution is that the determinant of the set 
(A. 23) equals to zero, that is 


h i SNX+h 4 CSX 


h^SN X +h 2 CSX 


SNA 


DILI 


n J CSX -hi SN X + CII X h$CSX -h 2 SN X -i SH X 


csx-cnx 


- SNX-SIIX 


‘AC II A 


_ jC . J/X + J^W(SIVX-SiJX) -2 SHX + ^^kl(CSX-CH X) | 

(A. 24) 
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Actually, Eq. (A. 24) is the characteristic equation of the Bemoulli-Euler beam equation used 
in this paper, from which the eigenvalues will be derived. Solving the characteristic equation 
by using binary searching technique, we can obtain the eigenvalues AT,-, 

K, — 0.0172, 0.0507, 0.1074, 0.1405, 0.1998, 


0.2258, 0.2263, 0.2268, 0.2278, 0.2283,... 

After determining the eigenvalues Kf s we can solve the set (A.23) for the coefficients 


C 3 , C±, di and d 4 . In fact, there are infinite number of solutions to the set (A.23) because 
C'i, C. u </ 3 and d 4 are not totally independent. Assuming the solution to the set (A.23) is 


' C 3 = a,C, 

Ca = PiCi 

< 

d$ — v |C| 

> d A = Ci 

then, from Eqs. (A.l 1) and (A. 18) we have the characteristic functions 


(A. 25) 


^ L t 00 — H (*hk t x — yin k t x) -f (chk t x — cos k t x) 

0<x ^ x i o 

5 1 ^ Y >i, i x ) - ( ‘« (/‘I (hi) sin k,.r f hi ( k,l ) cos k,l + shk,x j (A. 26) 

+A (k,i) sin k t x + li> ( k,l ) cos k,l + chk,x] 

^*10 5 ^ a : < / 


By superposition, then, the solution to the PDE should be 


i 

= * ‘ Mi tusw i ( + siiuu.ij (A. 27) 

1 

where, u>, - A'“ a. 

When proportional damping is taken into account, the PDE describing damped lateral 
vibration will be 


m 


d 2 y dy d A y 


Of- 


+ c ft + * 


a t * 


= 0 


(A. 28) 
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(A. 29) 


Here, two parameters a and b, which relate the coefficients of the PDB, are defined, 

< n 2 = i- = 

I ^ m pA 

l 26 = ^ 

The substitution of y (x,t) = (t)Vj(x) into Eq. (A. 28) yields 

i 

£ (mf.v- + C’T.r, + fcT.r"") = 0 

t 

which bears further a set of independent equations under the generalized coordinates T, (t) by 
the orthogonality property of the eigenfunctions, 

m,fi + c,T, + k t T t = 0 ( A - 3 °) 


where. 



Eq. (A. 30) can be expressed in modal form 


generalized mass 


generalized stiffness 
generalized damping 


where. 


f, + 2 &W U ,T, + U >lT x = 0 




c I 
m, 


(A. 31) 


Note that the damping ratio & is related to the eigenvalue K, through the parameter a and b 
In fact, 


Zi = 


1 £i 

2u; 1Ji m, 


1 c _ b 
2w„, m aL? 


(A. 32) 
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Eq. (A.31) turns out to be the equation of motion of a simple mass-spring-dashpot system. 
Therefore, the ith-component of the response can be expressed as 

Ti{t) — c (A, cos u>j t t t B, sinwj ( <) (A. 33) 

where, 


uj. 


(A. 34) 


ll ' u '"' Z 1 “ ^ ~ ~ = 'J ~ b ~ 

Thus, T,(t) can be expressed in terms of parameters a and b, 

T '(t) = e~ bl ^.4, cos t \J (uk'f) 2 - b 2 + Bi sin ty/(akf) 2 - 
By superposition, finally, the solution to the Eq. (A.28) should be 


i/(.M)- Y](.C )T,{t) 

l 

— ^ 1 , (.t) c bt ^4, cos \J (akj ) 2 - b 2 + B, sin tyj (aA.“) 2 - 6-^ (A. 35) 

where, K t (:r) are the eigenfunctions shown in Eq. (A. 26). 
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